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Abstract 

Assuming equation of state for quintessential matter: p = w(z)p, we analyse dynamical behaviour 
of the scale factor in FRW cosmologies. It is shown that its dynamics is formally equivalent to 
that of a classical particle under the action of ID potential V(a). It is shown that Hamiltonian 
method can be easily implemented to obtain a classification of all cosmological solutions in the 
phase space as well as in the configurational space. Examples taken from modern cosmology [14 



me 

3 



illustrate the effectiveness of the presented approach. Advantages of representing dynamics as a ID 
Hamiltonian flow, in the analysis of acceleration and horizon problems, are presented. The distant 
supernovae type la data are used to reconstruct the expansion scenario. The inverse problem of 
reconstructing the Hamiltonian dynamics (i.e. potential function) from the luminosity distance 
function cIl(z) for supernovae is also considered. 
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I. INTRODUCTION 



The Newtonian analogue to the Friedmann- Robert son- Walker models (hereafter FRW) 
was first considered by Milne and McCrea who used classical mechanics concepts to describe 
the expansion of the universe This approach has some difficulties connected with 

description of an infinite homogenous Newtonian cosmology Pressure, in contrast to 
Newtonian theory, has only relativistic character [jj. In this way, a quintessential universe 
has no Newtonian analogues. However, pressure seems to be play important role in the 
evolution of the real universe. Recently found accelerated expansion of our universe is due 
;o the presence of a certain vacuum energy with the negative pressure equation of state 



This type of energy is called dark energy and equation of state, in the general 
quintessential form, is usually postulated in the analysis of dynamics P|. 

Lima et al. [lOj considered the problem of reduction of FRW cosmologies with the 
equation of state p oc p to the problem of classical particle under 1-dimensional homogenous 
potential. We generalize this description for the FRW quintessential cosmologies. 

Recent observations of supernovae (SNIa) have revealed that this dark energy dominates 
the evolution of the universe. The most natural candidate to represent dark energy seems 
to be the cosmological constant. However, it is necessary to introduce a fine tunning of 120 
orders of magnitude in order to obtain the agreement with observations. Another popular 
idea is the so-called quintessence, a self-interacting scalar field jllj], but the quintessence 
program suffers from the fine tunning of microphysical parameters. In this work, we discuss 
the possibility that the dark energy is characterized by the equation of state 

p = w(a)p + = w{a(z))p + 0, (1) 

where z + 1 — a~ l is the relation between redshift z and the scale factor a(t) in FRW models. 
The main aim of this paper is to reduce the FRW dynamics (with the equation of state 



in the form ((TJ) to the form of the Newtonian equation o 



motion for a unit mass particle in 



ID potential. Therefore, we extend classical results j]]]-|4( to the quintessential cosmologies. 
On the other hand, the question concerning the reconstruction of the equation of motion 
(the same concerns reconstruction form of equation of state from the observational 
data becomes natural. We consider problem of reconstructing the potential function for the 
particle-like description of dynamics. 
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Organization of our paper is the following. In section 2, different methods of reducing 
the FRW quintessential cosmology (Q.C.) to the ID Hamiltonian flow is presented. Some 
applications of this formalism are given in section 3. In section 2, we present simple anal- 
ysis methods of classical cosmological problems such as the horizon problem, acceleration 
or cosmological constant problems. Finally, in section 4 the problem of reconstructing dy- 
namics (potential function) from (Il{z) in a model independent manner (as far as possible) 
is investigated. 



II. HAMILTONIAN DYNAMICS OF THE Q.C. MODELS 

We consider dark energy description in terms of single w(a(z)) but it is always possible to 
treat the equation of state of dark energy in terms of total pressure versus total density ratio. 
Moreover, in any case our treatment can be generalized to the case of many-component of 
non-interacting fluids. 

After assuming the quintessential fluid with the equation of state (JTJ), where p is the 
pressure and p the energy density, the dynamics of homogenous and isotropic models can 
be described by the following set of equations 

^ = -i(l + 3i»(o))p+p (2) 

£ = -3|& + «<«»* P) 
Equation (j2J) is the Raychaudhuri equation while eq. (jHJ) is the continuity equation for fluid 
with density p and equation of state ((TJ), ' = ^, t - cosmological time. 
The first integral of system ©-© is the Friedmann equation 

p = 3- + 3--A, 4 

where H = d(\aa)/dt is the Hubble function, A is the cosmological constant and k is the 
curvature index. 

After substitution (J3J) into (J2J) we obtain 

a + 'tp(a)d 2 + K,(a) = 0, (5) 



where 



1 + 3w(a) 

*P(a) = ^ , (6) 



(7) 

There are two different methods of reducting (J3J) to the form of the Newtonian equation of 
motion (i.e. of elimination of the ijj(a)a 2 type term from (J5J)). 
1. The reparametrization of the time variable 

t -> r : dt = (p(a(r))dr, (8) 



n(a) 



1 
2 



k 

(l + 3w(a))-- (l + w(a))Aa 



ip = exp J tp(a)da (9) 

reduces basic equation (jSJ) to 

„_d 2 a dV 
a =d^ = ~^ (10) 

where 

^(a) = J K(a)<p 2 (a)da (11) 

is the potential function for system (jlOj) . Finally, we obtain the dynamics reduced to the 
Hamiltonian flow in the 1-dimensional configurational space. 

H(a,d) = ^- + V(a). (12) 

It can be checked that 7i = E = const, is a constant of motion provided that energy density 
p satisfies continuity condition (jSJ). Therefore, trajectories of the system in the phase space 
(a, d) lie on the energy level 7i = E = const, (in the case of the vacuum cosmology we have 
E — 0). Of course, the Hamiltonian constraint should be consistent with the form of first 
integral (J3J). Hence, from (j3J) we obtain 

d 2 Tr . . d 2 k pa 2 Aa 2 

1- via) = 1 , 

2 w 2 2 6 6 

i.e. the potential for the system ()12j) is 

Now, the physical trajectories lie on the zero-energy level 7i = E = which coincides with 
the form of the first integral, because now the curvature and A-terms are included into 
the effective density p cS . Formally, the curvature term can be absorbed into the potential 



function by postulating the curvature fluid for which Wk = —1/3, pk = —3k/a 2 (as well as 
cosmological constant term, where p\ = —p\, /?a = A). 
From (|12j) and (|13|) we can see that both Hamiltonians 

W = ^|^| +V(q), H = ^\^ \ + V 2 (a)V(q) 




reproduce, in the consistent way the same Friedmann equations and they also reproduce 
equivalent equation of motion, because the gauge freedom in choosing the lapse function 
(i.e. freedom in reparametrization of time t — > r: dt = ipdr: r fi = y? 2 (a)7i). 
Let us consider special case of constant w. We have 

ip(a) = 
V(a) = 

\ z, u 
n^n = a 1+3w l \a 2 + V(a) 
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Classical equations of motion are 



4 = a=-A (W «), (14) 

where 

N = a 1+3w . 

The lapse function N plays the role of a Lagrange multiplier and upon its variation we 
obtain Hamiltonian constraint Tt = 0. 

The cases of constant quintessential coefficient w are important in applications. The 
forms of potential function for different kinds of matter and flat model are presented in 
Table HI Of course the presented formalism can be simply generalized to the case of any 
mixture of noninteracting multifluids with the equation of state Pi = WiPi Q|. For our 
purpose it is useful to put the dynamical equations into a new form by using dimensionless 
quantities 

*=-, T=\H \t, n h0 = ^ (15) 

<20 Pcr.O 

with 

a 2 

H — —, Per ft = 3iJ . 
(X 



TABLE I: The forms of potential function for different kinds of matter and flat model 





stiff mat. 


rad. 


dust 


string 


top. def. 


A 


phantom 


w 


1 


1/3 





-1/3 


-2/3 


-1 


-4/3 


V 


oc a~ 4 


oc a~ 2 


oc a^ 1 


const. 


oc a 


2 

oc a 


oc a 3 



The subscript means here that the quantity with this subscript is evaluated today (at time 
t ). Additionaly, we define = —"ik/QH^ and Q\ = A/3H$. Then the Hamiltonian is 

H=£ + V(z), (16) 

where 

V{x) = —^Qf.,o n ifi^~ 3wi (17) 

i 

which should be considered on the zero-energy level. 
The basic dynamical equations are then rewritten as 

y + V(x) = 0, 



x 

2 



i^a,o(l-3^)x- 3 ^. (18) 



In the general case, potential function can be obtained from p c g i.e. 
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where for the quintessence matter we have 



V{a) = - P -^±a\ (19) 



Pcfr = Pa + Pk + Poa J exp < - 3 / —^~da 

p cS = A - 3k{l + zf (20) 

- u,(a(*)) 



+ p (l + ^) 3 exp|3 J ^ 



+ z 



-dz >. 



2. In the second approach to the FRW dynamics representation by the Newtonian equation 
of motion, we redefine the position variable a, namely we define X in the following way 

a^X = a D(a \ (21) 



where D(a) is chosen in such a way that the term with d 2 is absent in (jSJ). Hence we obtain 
new variable 

/a r-a 
if (a) da, D(a) = \og a / if (a) da (22) 

and the equation of motion takes the form 

X = -K(a(X))if(a(X)) = (23) 
In the special case of constant w we have 

D(w) = log. ^a^da = log a a§( 1+ -) = 2(1 + w), 
and because of the obvious formula 

V(a(X)) = J K(a)if(a)dX, (24) 
in the considered case, we obtain 



f = a 1± ^, X = J^\ (25) 



H(X,X) = ^- + V(X)=0, 



= (i+„) x » - ? x 2 + ,„ 



(26) 



Let us observe that due to the new variable X, the considered system assumes the very 
simple form. It can be now considered on the energy level Ti, = E > 0, E — |(1 + w)po and 
V(X) = ~X 2 ( 1- d) — ^X 2 . This system was analyzed on the phase plane (X,X) many years 
ago from the point of its structural stability. It is interesting that there exists its natural 
generalization to the case of quintessential cosmology. 

Let us now emphasize some advantages of the considered approach. In general, the choice 

of 

X = P v^exp (- f 1 ^^da') da (27) 



makes it possible to reduce Q.C to the nonlinear particle mechanics. There are many 
advantages of using language of nonlinear mechanics. Let us consider some of them. 



1. Representation of dynamics as a one-dimensional Hamiltonian flow allows us to make the 
classification of possible evolution paths in the phase space as well as in the configuration 
space. Then the discussion of the existence and stability of critical points can be performed 
based on the geometry of the potential function. It is so because the stability of critical points 
is determined by the Hessian (d 2 T-C/dx l dy l ). In our case the Hamiltonian function assumes 
the very simple form characteristic for simple mechanical systems for which the Lagrangian 
function is natural, i.e. quadratic in velocities C = \g a pd a q 15 —V (q) . In our case, g a p = const. 
Then the characteristic equation for linearization of the system is A 2 + detA = 0, where Aj 
are eigenvalues of the linearization matrix A of the Hamiltonian system. Therefore, the 
only possible critical points in a finite domain of the phase space are centres (detA > 0) or 
saddles (detA < 0). Of course, we can explore dynamics given by the canonical equations. 
Denoting x = a, y = a = x we obtain 2D dynamical system in the form 

. &H 

x = -*- = y> 

dy 

. _ &H _dV 

^ dx dx 

We can observe that trajectories are integrable in quadratures. Namely, from the Hamil- 
tonian constraint 7i = E = we obtain the integral 

r da , . 

t-t = / ; 29 

and for some specific forms of the potential function we can obtain the exact solutions. 

2. It is possible to make the classification of qualitative evolution paths by analyzing the 
characteristic curve which represents the boundary equation in the configuration space. For 
this purpose we consider the equation of the zero velocity, a = 0, which represents the 
boundary of the domain admissible for motion T>e = {a G K+ : 1{E — V) ^ 0}. By 
considering the boundary of Ve given by the condition 

dV E = {a : V{a) = 0}, 

full qualitative classification of evolutional paths can be performed. 

3. We can find the domains of cosmic acceleration as well as the domains for which the 
horizon problem is solved. Because a = —dV/da, one can easily see that acceleration of 
the universe takes place if V(a) is a decreasing function of its argument. The condition 
dV / da\ a =a determines the static critical point on the phase plane. 
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Another interesting question concerns the horizon problem. It is easy to prove the 
following criterion of avoiding this problem. The FRW cosmological model does not have 
an event horizon near the singularity if a(t)c~ l tends to a constant while a(t) tends to zero 



When all events whose coordinates in the past are located beyond some distance du then 
can never communicate with the observer at the coordinate r = (in R-W metric). We can 
define the distance du as the past event horizon distance 



Of course, whenever / diverges as to ~~ * 0, there is no past event horizon in the spacetime 
geometry Q]. Then it is in principle possible to receive signals from sufficiently early 
universe from any comoving particle like a typical galaxy. It the t' integral converges for 
to — > then our communication with observer at r = is limited by what Rindler has called 
a particle horizon. The particle horizon will be present if energy density is growing faster 
then a~ 2 ~ e as a — > (e > 0) - Therefore if a > a~ e ^ 2 (d > A) then there is a particle 
horizon in the past. 

Let c~ l a < A (c = const, is assumed but the corresponding theorem can be established 
for the case of variable c(t) [3]). Then I ^ \ Jq° = ^(lna + oo). Therefore, / diverges as 
a — > and there exists no past horizon if the velocity of expansion factor a ^ A is bounded. 
Our investigation of the particle horizon is independent of any specific assumption about 
the behaviour of a(t) near the singularity or of specific form of the equation of state. If 
one assumes a linear equation of state p = wp and w = const., then Friedmann's equations 

2 

imply the following behaviour for a(t) ^ (t) 3(7+1) near the singularity t = 0. The integral 
I would thus diverge only if 7 < —1/3, i.e. only if the pressure p becomes negative. This 
is the condition for solving the horizon problem and it is identical to that for the solution 
of the flatness problem (8irGp/3 term will dominate the curvature term in a long time 
evolution). We can also show here that the integral J would diverge only if the pressure 
of the cosmological fluid takes negative values in the general case of w(a(t)). Conservation 
condition can be rewritten in the form 



121. 
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We can werify that boundedness of a(t) means also that pa 2 remains bounded near the 
singularity (see FRW eq.). Therefore, pa 3 — > as t — > 0. By integrating both sides of the 
equation written above from to t we obtain 

-3 / pa 2 adt = a 3 pit) ^ 0. 
Jo 

Consequently, w(a) must assume negative value without any specific assumption about the 
equation of state. 

The above criterion can be now formulated in the language of the phase space variables. 
If V(a) ^ const, (it does not diverge at singularity) as a — > then there exists no past 
horizon of the particle. 

4. Another natural generalization of the presented formalism consists in including anisotropy 
in simple Bianchi models (Bianchi I or Bianchi V models) or VSL models with variable 
velocity of light of the Albrecht-Maguejo-Barrow type. 

In the case of model BI (k = 0) or BV (k = —1), basic equation assumes the generalized 
form 

- = -l^- 1 -(l + 3w(a))p+^, (30) 
a 3 o 3 

where c(a) = Coa n (in the AMB parametrization) . 
Equation (|3"Uj) has generalized first integral 

2 



a \ p kc 2 a 2 Ac 



a 3 a 



2 



(31) 



where a 2 = \a a b& ab is the shear scalar which can be used to reduce dynamical equations to 
form (J3J), but now ip(a) preserves its form whereas n(a) — > R{a) takes the new one which 
builds the new potential function. 



c*(a) 



k 

'1 + 3w(a)) (1 + w(a))Aa 

a 



v 2 a,, (32) 



+ _(1_ U ,( G )) ! 



V(a) — > V(a) = j K,(a)(p (a) da. 

5. Let us now consider a class of FRW cosmologies with dissipation in the form of bulk 
viscosity. Owing to the assumed spacetime FRW symmetries the shear viscosity vanishes, 
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and we deal only with bulk viscosity dissipation. It was given by Weinberg who classified 
the physical significance of bulk viscosity effects within the framework of general relativity. 
The presence of bulk viscosity is equivalent to introducing effective pressure 

p cS = w(a(z))p - 3£H, (33) 

where £ is the bulk viscosity coefficient (constant for the sake of simplicity). 

If we introduce new variable X following previously described scheme, we obtain the 
generalization of (|23|) in the form 

dV 

X-aX + — = 0. (34) 

Equation (|34|) has the form of nonlinear oscillator equation with the damping force oc X 
and the exciting force F(X) = — In the following, this system will be analyzed qualita- 
tively by using methods of dynamical system. Therefore, on the basis of the corresponding 
reduction, the particle-like description may be given if effects of bulk viscosity are included. 
6. If we specialize the standard situation in cosmology with A = and fluid which satisfy 
weak energy condition p+p ^ in which we set p = wp and w = const, then the conservation 
equation gives p oc a~ 3 ^ 1+w ^ and the p/3 term will dominate the curvature therm ka~ 2 at 
large a so long as the matter stress violate strong energy condition p + 3p < and obey 
p + p ^ 0, that is if — 1 ^ w ^ —1/3. This is what we shall mean by the flatness problem. 

2 

The scale factor then evolves as a(t) oc t 3 <™+ 1 ) if w > — 1 or exp (H t), H constant if w = —1. 

In the general case of variable w(a(z)) quintessential model will provide a solution of 
flatness problem if zero curvature solution is representing late time attractor in the phase 
plane. This is possible in terms of potential function if d > const, as a » 1 or 

-C 2 /2 < V < (35) 

This is the condition oposite to that for the solution of horizon problem V < —C 2 /2 but it 
is formulate rather for large a than for small a. Let us notice that ifd>0asa>>l then 
a > Ct and period of accelerated expansion in long term behaviour can solve the flatness 
problem. 

Let us consider the case of non-zero cosmological constant added to the potential function. 
Then in order to explain why it does not totally dominate matter term p/3 at large a we 
would need the presence of fluid with subnegative pressure, i.e. with w < — 1 such that a 
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weak energy condition is violated p + p < 0. This is what we shall mean by the cosmological 
constant problem. If we assume that p is decreasing function of scale factor the solution of 
cosmological constant problem leads to contradiction a < 0. In the general case we have 
following condition for solving A-problem 

p cff >A^^(^) >0^\V\>Aa- 2 (36) 

i.e. modulus of potential function is growing faster than a~ 2 at large a. 

Let us notice that for w = const, above condition leads to the identical condition to 
that to solve the horizon problem — 1 ^ w ^ —1/3. It is consequence of the fact that 
corresponding conditions are valid at different domains a — > oo and a — > 0. 



III. APPLICATIONS 



We assume now that quintessential coefficient w(a(z)) = p/p is the analytical function 



of z 



Let 7q = w(0) and 



w{a(z)) = w(0) + ^ 
dz 



ld 2 w 



Z+ 2dz* 



z z + ... 



(37) 



Ti 



dw 

dz 



72 



ld 2 w 



2 dz 2 



i • • • i n 



li 



1 d l w 



i\ dz 1 



We consider models with vanishing A-term because conception of quintessence is treated as 
an alternative to the cosmological constant. The potential function of the dynamical system 
describing the FRW model with equation of state (|37|) is given by 



/ a k { f a 

K,(a)Lp 2 (a)da = — aexp < 3 / 

= -(l + z)- 1 oc P j(-3) J 



w\a) 



-da 



w(a(z)) 

Tz) 



(38) 



dz >. 



For example, in the case of mixture of noninteracting dust and fluid for which p x = w x p x 



where w x = const., we have 



,3w x 



V(a) = -a\ 



POm q3w x 



(39) 



12 



here we consider quintessential matter for which 



Px = Poxd 3 ( 1+Wx \ p 0x = const. 

and dust matter for which 

Pm = POmCt" 3 , pom. = Const. 

If for quintessential fluid w x = w x (a), then the potential function takes the more general 
form 



V(a) = \a 



exp (Sf^da) 



(40) 



^exp (^f^daj +1. ' 

After substituting (f3"Tj) into (jHHJ) we obtain the general form of the potential function 
which can be useful for a classification of cosmological models in the configurational space 



V(a(z)) = hl + z)- {l+ ^ 

f z ( 7l g + 72 ^ 2 + ■ ■ ■ ) 
x exp <( (-3) J ^—^ dz }. 



(41) 



Let us consider a few special cases of (}4~Tj) . 
A. 

Let 70 = — I (like for a string) which determines the boundary the strong energy condition 
violation. Additionaly, we consider expansion of w(z) up to the second order term. Then 
we obtain from the integral 

j llZ ^f dz = (71 - li)z + - (71 - 72) In (1 + z) (42) 

the exact form of the potential function 

V(a(z)) = - e -^~^>-h^\\ + z )3(7i-72) ; ( 4 3) 

or for any 70 



Of course, the trajectories of the considered system lie in the domain admissible for 
motion defined as 

V E = {aeR + : 2(E - V(a)) >0,E = 
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o / / wla) 
p = p$a ex P I 3 / da ] . 

Therefore, equation of the boundary dV E is 

y = ft <^ ^ e -3(7i-7 2 )2-§72^ 2 _ _j_ ^(l+37o)-3(7i-7 2 )_ 

If we substitute 70 = — |, then we obtain the potential function expressed in terms of a 
single parameter 71/72 = x, namely 



1 

1 + zY I 2E^ 



Finally, we obtain 



k 



\nE + z 2 /2 . , 

In (1 + zj — 2; 



where E = {2E/k) 1 ^ 2 . 
The plot of x (z) for different E is shown in Fig. |21 
Let us note that if po > (J5 > 0), then there is no boundary for k = —1, whereas if 
Po < (E < 0), there is no boundary for k = +1 of the domain admissible for motion. 
Finally, we consider the evolution path as a level of x — const, and then we classify all 
evolutions modulo the quantitative properties of their dynamics 
B. 

72 = 0, 70 7^ 0, i.e. we consider quintessential model with pressure p = 70P + (71 z)p. 
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FIG. 2: Diagram of dependence x(z) for classification of models with potential Q44|) 
(E = 0.1,0.6,1,2,10) 



Then we obtain the potential in the form (Fig 03) 



k 



V{a{z)) = -e~ 37l2 (l + ^)-( 1 + 3 7o)+37i. 



(46) 



Let us consider the boundary of the configuration space given by V(z). Then 7o can be 




FIG. 3: Diagram of the potential function V(z) for the case B, p = (70 + Jiz)p, 70 = —1/3, 
7l = -2/3/3 (J3 = -10,-1,0.1,0.4,1,10), k = ±1 



expressed as a function of z in the following way 



, , , 1 . v l\z c/3 

3 In 1 + z) In 1 + z) 



(47) 



where z = a 1 — 1, c = — In |^-|. The plot of 70 (z) for different 71 is shown in Fig. H]for po > 
and fc = +1. If we put 7o = —1 (cosmological constant) then we can obtain classification of 
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FIG. 4: Diagram of boundary curve 8T>e usefull in classification of the model in the configurational 
space (formula f37|>). We put here 71 = -2/3/9, {fi = -10, -5, -1, 1, 5, 10), c = -1 



cosmological models with cosmological constant in terms of levels of 71 = const., where 




FIG. 5: The plot of Ji(z(a)) (formula @$) for different signs of c (c = -5, -1.5, 0, 0.15, 2, 5, 10) 



c. 

p = w x p x , w x = const. Let us consider quintessential matter in the form of noninteracting 
mixture of dust and matter described by the equation of state p = w x p x , where w x = const. 
In this case, classification of all evolutional paths can be given in terms of potential function 
(Fig. El) 

k 1 a 3wx \ 

"W = 2« (3=5=+T ' (49) 

16 



where 

n — n n~ 3 ( 1+Wx ) n — n n~ 3 

Px — POxO j Pm — POmO- ■ 

The boundary curve V(a) = E can be used to classify all evolutional paths in the configu- 

±v 




0.5 1 1.5 2 2.5 3 



FIG. 6: Diagram of the potential function V(a) for the case C, p = w x p x , for different (w x = 

_4 _i _2 _l I U k _ 
3' - 1 ' 3' 3' x ' 3^' ^ — =C1 



rational space, namely 

In 



A 

a „ j4 P0m 



g (a) = 3w x (a) = — , (50) 

ma 



where A = 2E/k. 

The plot of g(a), for different A, is shown in Fig. [3 
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0.45/ 


V^O.55 
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FIG. 7: Diagram of g(a) for classification of evolution in the configurational space for the case of 
constant w x (A = 0.45,0.55,1,2), = 1 
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D. 

Since the observations of the supernovae of type la indicate that the universe must be 
today in an accelerated expansion the nature of the fluid responsible for such a behaviour 
has been object of many studies. While the most obvious candidate for such component is 
the vacuum energy the posibility that the dark energy might be described by the Chaplygin 
gas is seriously suggested [l^ . 

The Chaplygin gas has an interesting motivation connected with the string theory. If we 
consider a d-brane configuration in the d+2 Nambu-Goto action, the employment of the 
light-cone parametrization leads to the action of a Newtonian fluid with the equation of 
state p = —A/p, whose symmetries are the same as those of the Poincare' group. Hence, the 
relativistic character of the action is somehow hidden in the equation of state (for review see 
R. Jackiv, A particle field theorist's lectures on supersymmetric, non-abelian fluid mechanics 
and d-branes physics). 

Let us consider FRW model with fluid in the form of generalized Chaplygin gas for which 
equation of state is given by 

V = -j a , 0^a<l. (51) 

The energy-momentum conservation implies that Chaplygin gas density depends on the scale 
factor as 

/ Q \ 1+a 

P=[ A+ ^) J ' ^ 

where A and B are constants and the Chaplygin gas corresponds to the case a = 1; A is a 
positive constant because sound velocity of Chaplygin gas is v 2 s /c 2 = A 2 /p 2 . 

Equation (}52|) interpolates smoothly between a dust dominated phase, where p oc a~ 3 , 
and a De Sitter phase, where p = —p, through an intermediate regime described by the 
Zeldovich stiff matter p = p. 

Following our previous consideration, the dynamics is given by the hamiltonian 

H = y + V(a) = 0, 

where 

, « 1 „ 1 / „ B \ , k 
V(*) = - 6 P^ 2 = - 6 { A+ ^) « 2 +2- (53) 

The dependence of V(a) for a = 1 is illustrated on the Fig. |H] The domain admissible for 
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motion of the system with generalized Chaplygin gas is 

V = {a: V(a) ^ 0}. 

The boundary curve &Do can be used to classify possible evolution paths in the configu- 
rational space in the following way. We consider constant levels of A (a) relation given by 

., , f3k\ 1+ -" B 

A ( a )=[y) -JW (54) 
The zero velocity curve for k = +1 are shown on Fig. Eland Fig. ITUl 




FIG. 9: Diagram of A(a) for classification of possible evolution paths for the case of constant 
k = +1, q = 0.5 and different B (B = 0, 0.6, 1, 2, 5) 



It is interesting that the existence and character of the critical points of the considered 
system a = x, x = — ^ depends on the geometry of the potential function. 

19 



A 




FIG. 10: Diagram of A(a) for classification of possible evolution paths for the case of constant 
k = +1, B = 2 and different a (a = 0, 0.5, 0.75, 1) 




-3 
■4 

0.5 1 _ 15 2 2.5 3 

FIG. 11: The phase portrait (a, x) for case D (k = A = B = a = 1). The shaded region is a region 
of accelerated expansion of the universe. 




It can be easily to shown that in our case on a finite region of the phase plane (a, x) only 
saddle points are admitted because d 2 V/da 2 < at the critical point a = ao, dV/da\ a=ao = 0. 
Therefore, all points are hiperbolic (TrA = 0) and the system is structurally stable. 

In terms of V(a) the domains of accelerating trajectories can be easily found, namelly 
if dV/da < then the system starts to accelerate. Because the diagram of V(a) is upper 
convex, the static critical point will separate regions without acceleration from the domain 
in the phase space where trajectories accelerate. One can shown that the critical value of a 
is 
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At this point the diagram V(a) has the maximum. 



E. 

Interesting formulas for w(a) were already proposed by Corasaniti and Copeland Q|. 
They considered a broad class of tracking potentials for scalar fields, namely V(<f)) oc <ft~ a 
- inverse power low potential (INV) Q], V(<f)) oc <t>~ a exp (<j) 2 /2) - supergravity potential 



(SUGRA) (13, 1/(0) oc exp(-a0) + exp (/3<j>) (2EXP) (li 
Copeland- Nunes-Rosati model (CNR) [20]. 



The Skordis model 







and 



INV 
SUGRA - 
2EXP 
AS 
CNR 



-1.2 ' 

1e-05 



FIG. 12: The dependence of the equation of state factor w(a) versus the scale factor a/ao for 
different models of potential of scalar field. 



First, we aply our method to construct a potential function of the corresponding Hamil- 
tonian dynamical system. It assumes the following form 

■sr 1 \ Po _i ( q f 1 w(a) 
V (a) = a exp 3 / da 

6 \ Ja a 

= -£ exp j 3 (f, l^da + F 2 £ f -^da + F s (l - a) j j, (55) 

where 

/ ^ G = ^[ E i(-/?nx) r , m -Ei(l)], (56) 

^ a fl n=0 

where x = a — a c , (3 = 1/A, Ei - exponent integral function and coefficients F\, F2, F% are 
determined by the condition that w(a) takes on the respective values of w r , w m and w°, 
during radiation epoch (a = a r ), matter domination (a = a m ) as well as today (a — a — 1). 
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FIG. 13: The dependence of the potential of the Hamiltonian system against the scale factor o/ao 
for a different class of tracking potentials. 



The coresponding forms of w(a) function for a different class of potentials are illustrated 
in Fig. 

The function f r ,m{a) has the following form jlsT ] 



where the corresponding values of coefficients a r c ,rn are taken from Table I in ref |15| . 

After substitution formulas (JoTj) to (J55j) we obtain different forms of potentials for a 
different class of models (see Fig. H3|) . From this figure we can observe that in all cases a 
potential function is of the same type (upper convex), like for the Chaplygin gas. Therefore 
the phase portraits determined from the potential functions of the systems are topologically 
equivalent. From the physical point of view it means that considered models can be seriously 
treated as candidates for dark energy description. 

Let us notice that the dynamical system 




(57) 



dV 



a = x, 




can be transformed to (z,x) variables (a = (1 + z) *) 



and then we obtain 
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Therefore after the reparametrization of time along trajectories r — > 77 : drj = — (1 + z) 2 dr 

(now rj will be a decreasing function of time variable r) we obtain the dynamical system 

dz dx dV 

dr] ' dr] dz 

Of course, system (JH%|) can be analysed in terms of dynamical systems, i.e. in terms of the 
method of qualitative analysis of differential equations on the phase plane (z,x). System 
(j3Hj) has the first integral in the form 



^| -2(E-V(z)) 



0. 



(59) 



For the considered system only two types of critical points can appear, namelly centres or 
saddle points. If detA = ^t\ z _ zo is negative, the diagram of the potential function V(a(z)) 
has maxima; they correspond to the saddle point. On the other hand, if the diagram of 
the potential function V(a(z)) has minima, they correspond to centres. It is important 
that we can discuss the stability of critical points based only on the geometry of the po- 
tential function. It is easy to check that at the critical points of the system appearing at 

(|jj| = 0) we have k(o) = 0, and then 



a = const. (zq = const.), ^\ a 



d 2 V 



da 2 



2 3 k dw 
^ 2 a da 



3 ip 2 dw 
2 a da 



a=ao 



w(o)=-l/3 



3 9l .„ .ndw 



(60) 



z=z 



Therefore, if k^(z ) is positive, only saddles points can appear which guarantee the struc- 
tural stability of the system. 

For the case B ^ = 71, the above condition means that jik > 0. The phase portraits 



for 70 = —1/3, and various 71, are presented on Fig. IT9l 12012 fl and 1221 

Let us note that critical points are located at 

Zq = - 




(61) 



If we put, for example, 70 = —1 (cosmological constant term) then 71 > —| critical points 
can only exist on a finite domain of the phase plane because z > — 1. 

The idea of structural stability comes from Andronov and Pontriagin [3] . A dynamical 
systems S are said to be structurally stable if their dynamical behaviour remains qualita- 
tively (modulo homeomorphizm preserving orientation of trajectories) the same (equivalent) 
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under small perturbations. Structural stability is sometimes considered as a precondition 
for the "real existence". Structurally stable dynamical systems on the 2D compact space 
(for example on the plane with circle at infinity) form open and dense subsets in the space 
of all dynamical systems on the plane. Therefore, 2D structurally unstable models seem to 
be nonadequate for describing real proceses because of measurement errors. 

The main aim of the qualitative analysis of differential equations is not to find, and then 
to analyze, individual solutions but rather to investigate space of all possible solutions for all 
admissible initial conditions. A property is believed to be "realistic" if it can be attributed 
to large (rather typical than exceptional) subsets of models within the space of all possible 
solutions, or if it possesses a certain stability, i.e. if it is shared by a slighty perturbed model. 
There is a wide opinion among specialists that realistic models should be structurally stable, 
or even stronger, that everything that exists should possess a kind of structural stability. 

From the physical point of view it is interesting to answer the question: are the trajectories 
for which acceleration of the universe takes place, distributed in a typical or exceptional way? 
How are trajectories with interesting properties distributed in the phase plane? For example, 
the acceleration condition a = —dV/da > is satisfied if V(a) is a decreasing function of 
a. For us it is important that one should be easily able to deduce this from the geometry of 
the potential function only. 

In the phase space, the area of acceleration is determined by the condition that 

- = -i(l + 3uOp>0, ^ + ^(«)K) 2 <0. (62) 
a 6 oa 

(|62[) can be rewritten to the form which could be usefull in the analysis of the "probability 
of acceleration" which can be defined as a measure of the space of those initial conditions 
that lead to accelerating universes 

-Pn(E-V)+ 1 + 3w{a) <0, 
Oa a 

or 

g 

(1 + zf— In (E - V(z)) + (1 + 3w(a(z))(l + z) < 0. (63) 

For example, domains of acceleration in the configurational space, for the case D, as well 
as the domain of a > (QjL < 0) in the phase plane are presented on Fig. [TT] 
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Let us now consider the presence or absence of the particle horizon in the past. Good 
news from our earlier discussion is that this property can be detected from the shape of the 
diagram of the potential function of the system, namelly if V(z) goes to a constant (zero is 
included) as z — > oo then we obtain a model without the horizon. 

In the next section it will be demonstrated how we can answer the question about the 
horizon in the past on the base of V(z) taken from the observations. 



IV. INVERSE PROBLEM IN QUINTESSENTIAL COSMOLOGY 

The presented formalism gives us a natural base to discuss the redshift magnitude re- 
lation m(z) for SNIa supernovae observational data. But on the other hand, because the 
Hubble function is related to the luminosity distance, it is possible to determine both the 
quintessence parameter w(z) — p/p and the potential of the dynamical system V(z). There- 
fore, the equation of state as well as the whole dynamics can be reconstructed provided that 
the luminosity function di{z) is known from observations. It is called the inverse problem 
in dynamics of quintessential cosmology. 

As an example of constructing observables from the considered dynamics let us consider 
the luminosity-distance relation di(z) for quintessential models. 

If a light source of redshift z is located at a radial coordinate 7*1 (taken from R-W metric), 
its luminosity distance di, its angular diameter distance dA and its proper motion distance 
are given by 

d L {z) = (1 + z)a Q r 1 , d A = j^p~, d M = a ri, (64) 
where T\ calculated from metric gives 

da 1 f z dz' 



fin) 



(65) 



"o ad a Q J H(z') 

l + z 

sin -1 ri when k = +1 
7*1 when k = 

sinh 7*1 when k = —1 

Here ao is the present value of the radius of the universe. The above equation can also be 
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written in the form of a single compact equation as 



d L (z) 



l + z 




(66) 



where 



CO?) 
C(?) 
C(?) 



sing with «; = fi^o when fi^o > 0, 
sinh q with k = — f2&:,o when f^o < 0, 



with k — 1 



when f2 



fc,o — 



0. 



Thus for a given .M (absolute magnitude) and H(z) equation ()66|) gives the predicted value 
of m(z) (observed magnitude) at a given z. 

By using the K-corrected effective magnitudes which have also been corrected for the 
light curve width-luminosity relation and the galactic extinction, and by using the same 
standard errors <Ji,% and cr^f. of the supernova with redshift z% as used by Perlmutter et al. 
we compute \ 2 according to 



The best fit parameters are obtained by minimizing equation (|67jl. 

Luminosity distance and angular distance depend sensitively on the present densities of 
various energy components and their equations of state. For this reason, the magnitude- 
redshift relation for distant standard candles has been proposed as a potential test for 
cosmological models. 

In our formalism H(z) can be immediately taken from the first integral of the dynamical 
equation and then we obtain 



Formula (J68|) limits the determination of the luminosity distance because it depends on 
quintessential parameter w(z) through a multiple integral relation that smears out detailed 
information about w(z). If w(a(z)) can be expanded as the Taylor series following (|37|) then 
we obtain the simplest formula without the double integration 




(67) 




(68) 




(69) 
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Many authors [2l( assume that a quite accurate luminosity distance may be obtained and 
then examined to answer the question of whether the equation of state of the expanding 
universe can be determined uniquely. 

Our idea is more general. To determine the structure and evolution of an astrophysical 
system of the universe, the equation of state is usually necessary. By equation of state 
of the universe we mean the relation between the total energy density of cosmic matter 
and the total pressure. However, the equation of state relevant to the universe has not yet 
been established. Our idea is to reconstruct it from the form of the potential of the system 
V(a(z)). 

Let us consider, for simplicity, a flat model for which the Hubble parameter is related to 
luminosity distance by the relation 

-l 



H(z) 



d I dAz) 



(70) 



dz \ 1 + z 

Then it is possible to determine the quintessence parameter 

•M"i-1*<^>£(^)- < 71 » 

Here the term —1 is established from the condition that for stationary solution H(z) = const. 
-^2 ( ) w ih vanish but such a solution can appear on the phase plane (H,p) only as an 
intersection trajectory of the flat model and the boundary of the weak energy condition 
p + p = 0. 

Equation (fTTj) can be rewritten in the form 

p=-p-3£(z)H, (72) 

where 

2, r2 , , d 2 ( d L (z) 



^ z ) = -H\z){l + z) — 



3 y /K 'dz 2 \l + z 

The dependence p(H) manifests the presence of bulk viscosity effects in the model £(z) = 
— |J^=- Therefore, from di(z) we obtain the reconstruction of w(z) as a mixture of the 
cosmological constant term and bulk viscosity. However, in determining w(z) there is an 
inherent limitation because the value w(z) is poorly resolved and no usefull constraint can 
be obtained concerning its time variation. Of course, the value of w(0) and dw/dz with 
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accuracy 0.1 and 0.15 respectively will be obtained from SNAP3 but for the reconstruction 
of w(z) higher derivatives can be useful. If we measure ^r-(O), i ^ 2, then in principle, the 
reconstruction equation of the equation of state is possible, because then 

dw ld 2 w.. 9 

w(z) = w{0) + - {0)z+ -— {0)z > + ..., 



where 4-^(0) is determined from the recurrence formula, namelly 



d n w 
dz n 



(0) 



2 d r 



w(0) 



3dz r 
-14 



2 d n+1 
(lnm + -(l + z)f- 
v ; 3 V ' dz n + l 

o 

2 d 

3 dz 



(In if), 



(73) 



(In if). 



Then 10(2) can be obtained from H(z) in the following way 
H(z) 



m+l 

In iff 2) ^Vn^^ (In if) 



<i n u> 



(0) ^ w(z) 



El cf 10 . 
^(0K + M0). 



(74) 



At present this idea can not be realized but the idea of reconstructing dynamics requires the 
knowledge of V(z) which can be calculated from the Hamiltonian constraint 



V(a(z)) 



Pesa (z) 



6 



\H\z)a\z) 



-1 2 



(75) 



.(i + *)£(^)J 

On the other hand, the knowledge of V(z) gives us information how the horizon problem can 
be solved. Reconstruction of the potential function from SNIa data is presented on Fig ITS 
We obtain the plot of the potential function from fitting 

d L (z) 



l + z 



A + A X {1 + z') + A 2 {1 + z'Y + 



+ AM + z'f + AAl + z 



>\4 



(76) 



dz' 



function to the SNIa observational data. 

We can observe that as z — > 00, V(z) as calculated from (J75j) goes to —00. This can be 
treated as an empirical evidence of the presence horizon in the past. Let us note that other 
problems of the standard cosmology can be discussed analogously basing on information 
about V(z) and its geometry. 
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FIG. 14: The plot of the potential function V(z) has been obtained from equation (|75|) and from 
fitting formula (|76|) to the SNIa observational data 

V. CONCLUSIONS 

In this work a class of FRW models with quintessence matter is examined in the context 
of the present acceleration of the universe. Our results are the following. 

1. We have given a mathemathical background for discussing physical content of 
quintessential cosmology. Its dynamics is reduced to the dynamics of the unit mass 
point particle in ID potential. Then different physical properties, like acceleration of 
the universe, existence of horizon, can be formulated only in terms of the potential 
function of the system. The proof of the corresponding condition is quite general. In 
particular, it is independent of any specific assumption about the behaviour of the 
scale factor near the singularity (such as the assumption of power low behaviour) or a 
specific form of equation of state. 

2. The dynamics is formulated in the hamiltonian formalism and the full classification of 
possible evolutions in the phase plane as well as in configurational space is given. In 
the near future, it will be possible to obtain, from SNAP3, the exact value of 70 and 
71 appearing in the equation of state p = (70 + 712)/?, and then we could automatically 
answer the question about the horizon. 

3. The effectiveness of our treatment of dynamics of quintessential models in terms of 
single-particle mechanics is demonstrated for a broad class of tracking potentials. We 
obtain topological equivalence of the phase portraits (for this case) with the dynamical 
system obtained for the potential function reconstructed from SNIa data. 
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4. The idea of reconstructing dynamics instead of quintessential coefficient is considered. 
It is what we called inverse problem in quintessential cosmology. We demonstrate that 
the reconstructed potential function of the system produces the particle horizon in the 
past and solves the flatness problem. 
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FIG. 18: The phase portrait (z, x) for case A p = (70 + 71^ + / j2^ 2 )p (k = —1, 72 = 1) 




-3-2-10 1 2 

z 



FIG. 19: The phase portrait (z,x) for case Bp = (70 + Jiz)p (k = 1, 70 = —1/3, 71 = —2/3) 
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FIG. 20: The phase portrait (z,x) for case Bp = (70 + Jiz)p (k = 1, 70 = —1/3, 71 = 2/3) 
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